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ABSTRACT 


An  alternate  description  of  the  Plasma  Opening  Switch 
(POS)  is  presented.  In  addition,  a  preliminary  study  was 
conducted  into  the  ionization  of  desorbed  carbon  particles 
from  the  cathode  surface  in  the  switch  region.  The  formation 
of  a  locally  dense  plasma  layer  surrounding  the  cathode  is 
shown  to  influence  the  electric  field  and  the  Debye  screening 
distance.  This  effect  influences  the  emission  of  electrons 
from  the  cathode  as  well  as  providing  a  mechanism  for  opening 
the  switch. 
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A.  BACKGROUND 

Inductive  energy  storage  systems  are  capable  of  storing 
energy  at  densities  up  to  100  times  that  of  capacitive 
systems  [Ref.  1].  However,  inductive  storage  systems  are 
dependent  on  their  opening  switch  characteristics.  Switches 
range  in  type  from  mechanical  (slow)  to  exploding  fuses 
(faster)  to  plasmas  (fastest) . 

The  voltage  output  provided  by  an  inductor  is  given  by 
V  »  -L  di/dt,  where  L  [H]  is  the  inductance  and  di/dt  [A/s] 
is  the  changing  current .  As  the  switch  determines  the  di/dt 
term,  it  is  clear  that  faster  switch  performance  can  provide 
larger  voltage/power  output  with  smaller  input  current 
sources . 

Currently,  plasma  switches  are  being  incorporated  into 
vacuum  inductors  to  provide  the  fast  switching  necessary  for 
pulsed-power  applications.  Of  particular  interest  is  the 
Plasma  Opening  Switch  (POS) .  The  POS  has  been  used  in 
experiments  with  the  Gamble  I  [Ref.  1]  and  the  Gamble  II 
[Ref.  2]  accelerators  at  the  the  Naval  Research  Laboratory 
(NRL)  as  well  as  with  the  Pollux  and  Kalif  pulsed-power 
generators  [Ref.  3]  at  the  Kernf orschungszentrum  in 
Karlsruhe,  West  Germany.  The  results  of  these  experiments 
have  demonstrated  that  the  POS  was  capable  of  conducting  very 


large  currents  before  opening  in  less  than  10  ns.  Meger 
[Ref.  1]  reported  voltage  and  power  multiplication  of  2  and  4 
times  that  of  experiments  without  the  switch  while  Ottinger 
reported  similar  results  from  experiments  using  the  Gamble  II 
accelerator  [Ref.  2] .  The  POS  was  renamed  the  Plasma  Erosion 
Opening  Switch  (PEOS)  when  Ottinger  presented  his  theoretical 
model  of  the  switch  [Ref.  4] . 

The  Naval  Research  Laboratory  has  been  actively  using  the 
PEOS  to  provide  pulse  compression  for  their  experiments 
involving  a  plasma  radiation  source  (PRS) .  These 
experiments,  sponsored  by  the  Defense  Nuclear  Agency  (DNA) , 
were  undertaken  to  develop  a  source  which  simulates  the 
radiation  produced  by  nuclear  explosions.  The  results  of 
these  experiments  have  demonstrated  that  the  performance  of 
the  PRS  was  directly  affected  by  the  use  of  the  PEOS.  The 
PEOS  significantly  reduced  the  prepulse  and  compressed  the 
power  pulse  rise-time.  This  effectively  inhibited  the  growth 
of  instabilities  in  the  plasma  radiation  source  and  produced 
a  more  efficient  yield.  [Ref.  2] 

B.  THEORETICAL  MODEL  OF  THE  PEOS 

A  theoretical  model  of  the  PEOS  has  been  developed  by 
Ottinger  [Ref.  4]  consisting  of  four  phases:  conduction 
phase,  erosion  phase,  enhanced  erosion  phase  and  magnetic 
insulation  phase.  A  complete  description  is  available  in 
Reference  4. 
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During  the  May  1986  DNA  Pulsed-Power  Review  Meeting  in 
Las  Vegas,  Nevada,  the  Maxwell  Laboratory  presentation 
reported  the  results  of  their  latest  PEOS  experiments  on  the 
Blackjack'  5  generator.  The  results  demonstrated  PEOS 
switching  of  currents  up  to  2  mega  ampere  with  prepulse 
suppression.  However,  the  report  also  indicated  that  the 
switch  opened  too  slowly  for  pulsed  power  applications  and 
that  the  magnetic  insulation  of  the  cathode  surface  was 
incomplete.  The  conclusions  stated  that  critical  aspects  of 
the  switch  parameters  and  the  physical  mechanisms  of  the  PEOS 
model  were  not  well  understood.  [Ref.  5] 

Schwirzke  has  asserted  that  the  ion  diode  modeling  of  the 
PEOS  is  not  applicable  and  the  physics  describing  the  plasma 
switch  are  incomplete.  He  has  proposed  an  alternate 
description  of  the  physical  processes  of  the  switch  based  on 
reduced  field  emission  of  electrons  from  the  cathode  surface. 

C.  PLASMA  SHEATH  MODEL 

The  proposed  model  is  based  on  Schwirzke 's  earlier  work 
[Ref.  6]  and  the  recent  work  of  Halbritter  [Ref.  7] .  These 
experiments  demonstrated: 

*  Plasmas  in  contact  with  metals  induced  significant 
desorption  of  gases  and  other  impurities  from  the 
surface. 

*  Electron  field  emission  at  contaminated  surfaces  is 
achieved  with  reduced  electric  fields. 

These  observed  results  led  Schwirzke  to  postulate  an 

alternate  description  to  the  PEOS  model. 
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Schwirzke  has  asserted  that  Ottinger's  PEOS  model  does 
not  fully  describe  the  physics  of  the  switch.  Schwirzke' s 
major  points  of  disagreement  are: 

*  The  electron  current  is  field  emitted  rather  than  space- 
charge-limited  . 

*  The  sheath  that  forms  between  the  bulk  plasma  and  the 
cathode  surface  adjusts  to  the  applied  voltage  at  the 
cathode  rather  than  being  fixed. 

*  There  is  no  physical  mechanism  which  causes  the  PEOS  gap 
to  increase. 

D.  PURPOSE  OF  THE  RESEARCH 

The  purpose  of  this  research  is  to  qualitatively  describe 
Schwirzke' s  opening  switch  model  and  perform  some  preliminary 
testing  to  explore  the  influence  of  the  neutral  carbon 
ionization  on  the  performance  of  the  switch.  Other  aspects 
of  the  model  will  be  treated  qualitatively.  It  is  hoped  that 
this  research  will  serve  as  a  basis  for  the  continued 
development  of  this  model  and  that  important  insight  into  the 
switching  behavior  of  the  plasma  can  be  achieved. 


II.  PLASMA  SHEATH  i 


A .  OVERVIEW 

1 .  Purpose  of  the  Model 

The  proposed  plasma  sheath  model  is  an  alternate 
description  of  the  switch  model  known  as  the  Plasma  Erosion 
Opening  Switch  which  has  been  developed  by  Ottinger  [Ref.  4] . 
The  plasma  sheath  model  is  based  on  the  reduced  field 
emission  of  electrons  from  a  contaminated  cathode  surface. 
This  model  attempts  to  simplify  and  gain  insight  to  the 
understanding  of  the  physical  behavior  of  a  highly  non- 
uniform  plasma  within  a  confined  switch  region. 

2 .  Opening  Switch  Description 

The  basic  configuration  of  experiments  using  the 
plasma  opening  switch  (POS)  is  shown  in  Figure  1.  It 
consists  of  a  generator  section,  vacuum  coaxial  inductor 
section,  switch  section  and  load  section.  Plasma  is  injected 
from  plasma  guns  located  off-axis  toward  the  cathode.  After 
a  delay  period,  the  generator  is  fired  and  current  begins 
flowing  from  the  outer  anode  to  the  inner  cathode  through  the 
plasma  switch  region.  The  switch  region  conducts  the  current 
flow  for  some  period  of  time  before  opening  quickly  and 
providing  a  high  voltage  pulse  to  the  load  section.  It  has 
been  shown  that  the  performance  of  the  opening  switch  is 
highly  dependent  on  the  plasma  density,  the  delay  time  for 
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firing  the  generator  and  the  physical  dimensions  of  the 
switch  region.  [Ref.  8] 


Figure  1.  Experimental  Setup 


B.  PLASMA  SHEATH  DESCRIPTION 
1 .  Assumptions 

The  major  assumptions  which  apply  to  this  model  are: 

*  The  anode-cathode  potential  difference  is  applied  over 
the  sheath  region.  The  effective  resistance  of  the  bulk 
plasma  is  very  small  and  it  is  determined  by  the  Spitzer 
resistivity  and  the  switch  region  dimensions.  The 
potential  difference  required  to  drive  the  current 
through  the  bulk  plasma  is  minimal  and  can  be  neglected 
compared  to  the  total  anode-cathode  potential 
difference. 

*  The  stainless  steel  cathode  has  several  monolayers  of 
contaminants  on  the  surface.  The  main  source  of  the 
contaminants  is  from  the  carbon  plasma  guns.  After  each 
firing  of  the  plasma  guns,  a  certain  amount  of  residue 
settles  on  the  surface  of  the  cathode.  For  the  purposes 
of  this  research,  only  the  carbon  neutrals  will  be 
considered. 


*  The  streaming  gun  plasma  induces  desorption  of  the 
cathode  surface  contaminants.  The  desorbed  contaminants 
expand,  forming  a  dense  layer  or  cloud  of  neutral 
particles.  These  particles  are  ionized  by  the  emitted 
electrons  from  the  cathode  surface. 

*  The  initial  current  flow  in  the  switch  region  is  uniform 
and  in  the  radial  direction.  The  field  emitted 
electrons,  normally  emitted  at  very  small  sites  such  as 
whiskers  or  dirt  spots,  are  considered  to  be  distributed 
over  the  entire  cathode  surface  in  the  switching  region. 
In  addition,  the  ion  flow  is  considered  to  be  uniform 
with  a  constant  drift  velocity  toward  the  cathode. 

*  The  switch  dimensions  and  the  density  of  the  bulk  plasma 
remain  constant  throughout  the  duration  of  the  switching 
period.  The  density  of  the  bulk  plasma  is  confined 
without  losses  through  any  boundaries  with  the  exception 
of  the  ion  current  at  the  cathode  surface.  However, 
this  loss  is  small  and  can  be  neglected  so  that  the 
plasma  density  is  essentially  constant. 

*  The  sheath  region  does  not  break  down  due  to  the 
voltages  across  the  gap.  In  addition,  the  plasma 
temperature  assumes  the  energy  of  the  cathode.  As  the 
cathode  energy  increases,  the  sheath  distance  increases 
by  pushing  the  electrons  further  away.  This  imparts 
energy  to  these  electrons  which  in  turn  is  applied  to 
the  bulk  plasma. 

*  The  cathode  is  cold  so  that  thermal  emission  of 
electrons  is  not  significant  and  can  therefore  be 
neglected. 

2.  Desorption 

As  the  injected  plasma  comes  into  contact  with  the 
cathode,  the  neutral  contaminant  layers  are  disrupted  and 
these  atoms  are  desorbed  from  the  surface.  These  neutrals 
expand  slowly  (approximately  330  m/s)  and  form  a  locally 
dense  layer  of  atoms  to  be  ionized.  The  expansion  of  the 
neutral  particles  is  considered  only  during  the  time  period 
that  the  plasma  is  in  contact  with  the  cathode  surface  prior 
to  the  firing  of  the  generator.  This  is  due  to  the  slow 


expansion  velocity  and  very  short  time  frame  for  current 
conduction.  Therefore,  the  number  of  neutral  particles 
desorbed  and  the  distance  traveled  from  the  cathode  surface 
is  considered  constant  throughout  the  switching  phase. 

3.  Initial  Field  Emission 

The  initial  emission  of  electrons  from  the  cathode 
surface  is  produced  by  localized  spots  with  high  electric 
field  strengths.  The  electric  field,  E  (V/m)  is  enhanced  by 
a  factor  15  such  that  the  enhanced  field  is  given  by  F  =  J5E. 
The  enhancement  factor  range  is  10  £  15  £  1000. 

Halbritter's  work  [Ref.  7]  demonstrated  that  the 
electron  field  emission  current  from  contaminated  surfaces 
was  observed  at  electric  field  strengths  on  the  order  of 
F  =  107  V/m.  The  contributing  factors  to  field  enhancement 
are  dust  particles,  whiskers  and  other  contaminants. 

The  electric  field  is  defined  as  E  =  -  V<>,  where  <|>  [V] 

is  the  potential  field.  For  the  case  of  cylindrical 
coordinates  and  a  potential  field  which  only  varies  in  the 
radial  direction,  then  E  =  -d<j>/dr.  Furthermore,  if  it  is 
assumed  that  d<|>/dr  is  constant  within  a  given  region,  then  E 
can  be  approximated  by  the  potential  difference  over  that 
region . 

The  plasma  sheath  model  assumes  that  the  anode-cathode 
potential  difference  occurs  within  a  sheath  region  given  by 
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eoe<fc 

V  “ r~  (i) 

ne 

where  eo  [F/m]  is  the  permittivity  of  free  space,  e  [C]  is  the 
elementary  charge,  <(>c  [V]  is  the  cathode  potential  and  n  [rtf3] 
is  the  electron  density  of  the  plasma.  This  is  the  Debye 
screening  distance  for  a  negatively  biased  cathode.  As 
stated  earlier,  E  ~  <j)/XD  so  that  E  varies  with  the  plasma 

density  and  cathode  potential  as 


4 .  Ionization 

The  ionization  rate  of  the  neutral  particles  is 
dependent  on  the  neutral  density,  ionization  cross  section 
and  the  number  of  incident  electrons.  This  research 
considers  only  the  neutral  carbon  atoms  which  are  ejected 
from  the  surface  for  calculations  of  the  ionization  cross 
section . 

As  the  neutrals  become  ionized,  a  locally  more  dense 
plasma  forms  near  the  cathode  surface  which  influences  the 
electric  field  and  the  sheath  width.  From  Equation  2,  it  is 
clear  that  for  a  constant  potential,  the  electric  field 
varies  with  the  density.  For  increasing  density,  the 
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electric  field  strength  increases  and  additional  field 
emission  from  the  cathode  is  possible. 

5.  Ion  Removal 

The  ionization  of  the  neutral  particles  produces  a 
locally  more  dense  plasma  near  the  cathode  surface.  The 
ionized  particles  in  the  sheath,  under  the  influence  of  an 
enhanced  electric  field,  are  accelerated  toward  the  cathode 
surface.  The  time  required  for  these  ions  to  reach  the 
surface  is  determined  by  the  accelerating  field  and  the 
distance  traveled.  As  the  carbon  ions  are  removed  from  the 
sheath  region,  they  contribute  to  the  total  ion  current  flow 
due  to  the  higher  density.  In  addition,  the  reduction  in  the 
local  charge  density  causes  the  sheath  width  to  increase  and 
the  electric  field  to  decrease.  As  the  magnetic  field 
strength  increases,  the  electrons  begin  to  drift  toward  the 
load  side  of  the  switch  before  crossing  the  gap  and  entering 
the  bulk  plasma. 

6.  Magnetic  Insulation 

Magnetic  insulation  occurs  when  the  magnetic  field 
near  the  cathode  is  sufficiently  large  to  cause  the  electron 
Larmor  radius  to  be  smaller  than  the  sheath  width.  As  the 
magnetic  field  experienced  by  the  emitted  electrons 
increases,  the  electrons  experience  the  ExB  (E  is  assumed  to 
be  radial  and  B  is  in  the  0  direction)  drift  which  causes  the 

electrons  to  drift  parallel  (z  direction)  to  the  cathode 
surface  for  the  length  of  the  switch.  The  current  is 


switched  to  the  load  and  the  effective  switch  resistance 
increases.  The  switch  is  considered  open  when  the  Larmor 
radius  becomes  less  than  the  sheath  ”idth. 
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HI.  IONIZATION  OF  NEUTRAL  PARTICLES 


There  are  four  critical  areas  of  interest  in  the  plasma 
sheath  model  which  need  to  be  addressed.  These  are: 

*  The  ionization  of  neutral  particles  which  are  desorbed 
from  the  cathode  surface. 

*  Electric  field  strength  and  sheath  width  response  due  to 
the  increased  ionization  of  neutral  particles. 

*  Ion  removal  occurring  at  the  cathode  surface. 

*  Magnetic  insulation  which  causes  the  switch  to  open. 

This  section  is  devoted  to  the  calculations  necessary  to 

determine  the  number  and  density  profiles  of  the  desorbed 
carbon  atoms,  the  number  and  density  profile  of  the  newly 
formed  plasma  and  the  electric  field  produced  by  the  pressure 
force  associated  with  a  density  gradient.  The  computations 
used  are  based  on  fundamental  plasma  physics  relationships. 

No  effort  has  been  made  to  apply  statistical  methods  to  this 
problem. 

A.  DESORBED  NEUTRALS 

The  cathode  used  in  the  experiments  at  the  Naval  Research 
Laboratory  and  at  the  Kernforschungszentrum  in  Karlsruhe  is 
made  of  stainless  steel.  Most  stainless  steels  consist  of  a 
small  percentage  of  nickel  and  chromium  (usually  below  20%) 
mixed  in  with  the  base  metal  of  iron.  At  temperatures  below 
912  degrees  Celsius,  the  crystal  lattice  structure  of  iron  is 
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body  centered  cubic  (BCC)  with  a  lattice  parameter  of 
2.866xlO~10  meters.  The  contaminants,  which  are  situated  on 
the  surface,  will  occupy  the  interstitial  sites.  Therefore, 
the  BCC  configuration  and  lattice  parameter  of  iron  will  be 
used  for  the  calculation  of  the  number  of  carbon  particles 
desorbed  from  the  cathode  surface.  [Ref.  8] 

The  volume  of  a  layer  on  the  cathode  surface  within  the 
switch  region  is  given  by 


V  =  2nR  1  *  dr  (3) 

where  Rc  im]  is  the  radius  of  the  cathode,  1  [m]  is  the  length 
of  the  switch  region  and  dr  [m]  is  the  depth  of  the  surface 
layer.  The  value  of  dr  is  given  by  dr«0 . 5*m*2 . 866xlO-10, 
where  m  is  the  number  of  monolayers  desorbed.  The  BCC 
lattice  contains  two  atoms/cell,  therefore  the  number  of 
atoms  desorbed  is  given  by 

2nR  lm 

n= - S“in  (4) 

( 2.866x10  ) 


B.  NEUTRAL  DENSITY  PROFILE 

The  radial  expansion  of  the  neutral  carbon  layer  is 
determined  by 


(5) 
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where  Rg  [m]  is  the  off  axis  location  of  the  plasma  guns,  t 

[s]  is  the  delay  time  from  plasma  gun  firing  to  generator 
firing,  vd  [m/s]  is  the  plasma  drift  velocity  and  vn  [m/s]  is 

the  neutral  expansion  velocity.  The  neutral  particle 
expansion  distance  will  be  considered  constant  for  the 
duration  of  the  switching  action  due  to  the  short  time  frame 
involved  and  the  relatively  slow  expansion  velocity  of  the 
neutrals . 

Sub-regions  are  formed  by  dividing  Equation  5  by  the 
number  of  sub-regions  desired.  Initially,  each  sub-region  is 
assigned  an  equal  number  of  particles  based  on  the 
calculation  of  Equations  4  and  5.  The  profile  can  be 
adjusted  so  as  to  weight  the  number  of  neutrals  in  a  given 
sub-region  as  will  be  explained  in  Appendix  A.  From  this 
number  profile,  a  neutral  density  profile  for  each  sub-region 
will  be  calculated  from 


N„ 


*1<R*-R?) 


(6) 


[ 

> 
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where  NR  [m  x]  is  the  number  of  particles  per  sub-region  and 
R0  [m]  and  Rx  [m]  are  the  outer  and  inner  radii  of  the  given 
sub-region.  This  neutral  density  profile  is  important  for 
determining  the  mean-f ree-path  of  the  electrons  within  each 
sub-region . 
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C.  ELECTRON  FLUX 


The  number  of  electrons  flowing  through  the  layer  of 
desorbed  neutrals  is  determined  from  the  electron  current 
density  flowing  through  the  cathode  surface.  For  a  given 
total  current,  I  [A],  and  a  given  ion  current, 


■r  ",CVA 


(7) 


where  Ac  [m2]  is  the  cathode  surface  area,  the  electron  flux 
at  the  anode  is  determined  by 


*  10'9  (8) 

where  Aa  [m2]  is  the  anode  surface  area  and  T  has  the  units 
[number  of  electrons  per  square  meter  per  nanosecond] .  The 
electron  current  is  extracted  from  the  anode  surface  while 
the  ion  current  is  extracted  at  the  cathode  surface. 

Conservation  of  charge  requires  that  the  electron  current 
extracted  at  the  anode  must  equal  the  electron  current  at  the 
cathode.  Therefore,  the  radially  outward  electron  flux  into 
a  cylindrical  region  is  given  by 


r.  =  r 

tn 


(9) 
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where  Ri  [m]  is  the  radius  of  the  inner  surface  that  the 
electrons  flow  through.  The  number  of  incident  electrons  per 
nanosecond  into  each  cylindrical  region  can  be  expressed  by 


(10) 


which  is  constant  for  a  given  current. 

D.  DENSITY  OF  IONIZED  PARTICLES 

The  mean-f ree-path  of  the  electrons  is  determined  by 


X  =  (nno)'1  (11) 

where  <J  [m2]  is  the  electron  impact  ionization  cross  section 
for  the  given  target  material.  The  value  of  o  is  dependent 
only  on  the  energy  of  the  incident  electrons .  Table  1  is 
used  to  determine  the  electron  impact  ionization  cross 
section  for  neutral  carbon  atoms  [Ref.  9] .  The  values  used 
range  from  relatively  low  energies  ~  20  eV  up  to  ~  1000  eV. 

It  is  assumed  that  the  ionization  occurs  at  low  voltages  so 
only  relatively  low  electron  energies  will  be  considered. 

The  values  of  energy  are  given  in  electron  volts  and  the 
cross  section  is  given  in  square  meters. 

The  ionization  of  the  neutral  particle  cloud  is  dependent 
on  the  electron  flux,  mean-f ree-path  and  the  thickness  of  the 
medium.  For  a  given  flux  of  electrons  incident  upon  a  medium 


of  thickness  x  and  mean-f ree-path,  the  flux  of  electrons 
exiting  the  medium  is  reduced  by  the  factor  exp  (-x/X) .  The 

difference  between  the  incident  and  exiting  electrons 
represents  the  ionization  that  occurs  in  the  material. 


TABLE  1 

ELECTRON  IMPACT  IONIZATION  CROSS  SECTION 


Energy 

Cross  section 

Energy 

Cross  section 

20 

1.00E-20 

200 

1.50E-20 

30 

1.80E-20 

300 

1.20E-20 

40 

2 . 10E-20 

400 

9.00E-21 

50 

2.30E-20 

500 

8.00E-21 

60 

2 . 40E-20 

600 

7.00E-21 

70 

2.30E-20 

700 

6.50E-21 

80 

2.20E-20 

800 

6.00E-21 

90 

2.15E-20 

900 

5.50E-21 

100 

2.10E-20 

1000 

5.00E-21 

The  number  of  neutral  carbon  atoms  that  are  ionized  can 
be  expressed  by 


N.=  Ne[  1  -exp(-x /X)] 


(12) 


It  follows  that  the  density  of  the  desorbed  carbon  atoms  that 
become  ionized  is  given  by 


N. 


n- 

1 


*1  (Ro-  Ri> 


(13) 


At  the  end  of  each  nanosecond  time  increment,  the  number 


of  particles  that  become  ionized  in  each  region  is  subtracted 
from  the  original  number.  The  new  densities,  mean-free- 


paths,  etc.  are  then  recalculated  for  the  next  nanosecond 
increment . 

E.  DENSITY  RATIO  AND  PRESSURE  GRADIENT 

The  ratio  of  the  ionized  carbon  particles  to  the  plasma 
density  is  given  by 


d  =  T  <14> 

where  nt  [m-3]  =  nt  +  n,  is  the  total  ionized  density.  The 
electric  field  at  the  cathode  and  the  sheath  width  are 
influenced  by  increased  values  of  D.  The  electric  field 
experienced  at  the  cathode  can  be  modified  by  multiplying 
Equation  2  by  D1/2  which  demonstrates  the  influence  of  the 
ionization  on  the  electric  field. 

In  addition  to  the  influence  of  D  on  the  electric  field 
experienced  at  the  cathode,  an  increasing  D  produces  a 
pressure  force  within  the  plasma.  Treating  the  plasma  as  an 
ideal  gas,  the  pressure  can  be  expressed  as  p  =  nkT,  where  k 
[J/K]  is  the  Boltzmann  factor  and  T  [K]  is  the  electron 
temperature  of  the  plasma.  Considering  only  the  pressure 
gradient  force,  the  electric  field  produced  by  the  gradient 
can  be  expressed  by 

kT 

E  =  -  — —  Vn  (15) 

end 
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where  nd  [m-3]  is  a  modified  density  term.  This  term  is  the 
average  density  between  adjacent  regions  and  is  given  by 


n  +  n 
n.=  — — 
d  2 


(16) 


If  the  pressure  gradient  force  is  sufficiently  strong,  the 
motion  of  the  ions  and  electrons  in  this  region  will  be 
affected . 
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The  results  of  this  study  are  presented  in  two  parts. 
First,  the  ionization  occurring  in  the  desorbed  carbon  cloud 
is  evaluated.  Second,  the  effect  of  the  pressure  gradient  on 
the  electric  field  due  to  the  ionization  is  considered.  The 
switch  parameters  and  the  input  current  used  in  this  study 
are  taken  from  the  Gamble  I  experiments  reported  in  Reference 
10.  The  effect  of  varying  the  switch  parameters  on  the 
ionization  and  resultant  electric  field  strength  is  also 
evaluated . 

Typical  switch  parameters  given  in  References  1,  2,  4 
and  10  are  shown  in  Table  2. 

TABLE  2 

SWITCH  PARAMETERS 


Gun  Plasma  Density 

3 .  OxlO19 

m-3 

Plasma  Drift  Velocity 

7 . OxlO4 

m/s 

Delay  time 

1.5xl0'6 

s 

Time  of  flight 

1.36x10" 

6  s 

Cathode  radius 

2 . 5xl0"2 

m 

Anode  radius 

5.  OxlO"2 

m 

Gun-cathode  distance 

9 . 5xl0"2 

m 

Switch  length 

4  .OxlO"2 

m 

A.  IONIZATION  RATIO 

The  switch  parameters  listed  in  Table  2  which  will  be 
varied  during  the  evaluation  are  the  density  of  the  plasma 
guns  and  the  delay  time  between  the  firing  of  the  plasma  guns 


and  the  firing  of  the  generator.  In  addition,  the  effect  of 
the  electron  energy,  monolayers  desorbed  and  the  density 
profile  of  the  desorbed  particles  will  be  considered. 

1.  Density  profile 

The  density  profile  utilized  is  explained  in 
Appendix  A.  Essentially,  the  expanding  carbon  cloud  is 
divided  into  5  concentric  annular  regions  of  equal  spacing. 
Each  region  is  assigned  an  equal  number  of  neutral  carbon 
particles.  This  method  produces  a  slightly  non-uniform 
density  profile  due  to  the  cylindrical  geometry.  However, 
the  profile  can  be  altered  by  applying  a  weighting  factor  f. 
For  f  =  1,  the  number  of  particles  in  each  region  is  equal. 
For  0.5  <  f  <  1,  the  number  of  particles  is  more  heavily 
weighted  to  the  outer  regions  while  1  <  f  <  2,  more  particles 
are  allocated  within  the  inner  regions. 

The  effect  of  this  particle  weighting  can  be  seen 
in  Figures  2,  3,  4  and  5.  For  50  eV  electrons,  the 
ionization  of  the  carbon  particles  saturates  within  40-45 
nanoseconds.  Allocating  more  particles  in  the  outer  regions 
produces  a  less  varied  plasma  density  difference  between 
regions  than  allocating  them  to  the  inner  regions.  By 
weighting  the  inner  region,  an  order  of  magnitude  plasma 
density  difference  can  be  achieved. 

2.  Desorbed  layers 

As  seen  in  Figure  6,  variations  in  the  number  or 
fractional  number  of  monolayers  desorbed  affects  the 


ionization  ratio  by  the  same  factor.  Increase  the  number  of 
monolayers  desorbed  by  a  factor  of  10,  and  a  factor  of  10 
increase  in  the  ionization  ratio  is  observed.  This 
information  provides  a  baseline  to  consider  for  the  number  of 
monolayers  that  must  be  removed  in  order  to  produce  any 
significant  effect  on  the  electric  field  or  the  sheath  width 
as  proposed  by  Schwirzke.  Equations  1  and  2  scale  with  the 
square  root  of  the  plasma  density.  Figure  6  clearly 
demonstrates  that  in  order  to  influence  the  electric  field, 
at  least  0.1  monolayers  of  carbon  particles  must  be  desorbed 
from  the  surface.  In  order  to  see  the  trends  more  clearly, 
one  monolayer  will  be  utilized  for  all  future  comparisons. 

3.  Electron  energy  variations 

From  Table  1,  it  is  clear  that  electrons  with 
energies  below  100  eV  have  larger  ionization  cross  sections 
than  those  of  higher  energies.  Figure  7  demonstrates  that 
the  lower  energy  electrons  ionize  the  carbon  layer  much  more 
rapidly  than  the  higher  energy  electrons.  Based  on  the  known 
time  scales  of  the  plasma  opening  switch,  it  can  also  be 
inferred  that  the  voltage  across  the  sheath  must  remain 
relatively  low  (  20  £  f  ^  100  )  during  the  initial  conduction 
phase  of  the  current.  This  allows  for  rapid  plasma  density 
increase  near  the  cathode  surface.  However,  the  electric 
field  in  Equation  2  is  dependent  on  the  density  and  the 
applied  voltage.  If  the  anode-cathode  potential  difference 
remains  low,  then  the  density  must  increase  significantly  to 


provide  the  electric  fields  necessary  to  sustain  electron 
emission. 

4.  Delay  time 

By  varying  the  delay  time  between  the  firing  of  the 
plasma  guns  and  the  generator,  the  size  and  density  of  the 
desorbed  carbon  cloud  are  modified.  This  results  in  an 
overall  decrease  in  the  ability  of  the  electrons  to  ionize 
the  carbon  atoms  as  the  delay  time  increases.  As  stated 
previously,  the  effect  of  the  neutral  carbon  density  on  the 
mean  free  path  is  significant.  Figure  8  demonstrates  this  as 
the  density  ratio  D  decreases  as  the  delay  time  increases. 

5.  Plasma  gun  density 

As  the  plasma  gun  density  varies,  the  ionization 
ratio  also  varies.  Increased  plasma  gun  densities  require 
that  the  ionization  of  the  carbon  atoms  be  much  greater 
before  the  effect  is  noticeable.  Whereas  decreased  plasma 
gun  densities  can  produce  a  significant  ionization  ratio 
almost  immediately.  Figure  9  demonstrates  this  principle. 

B.  PLASMA  DENSITY  GRADIENT  EFFECT 

The  plasma  density  gradient  produces  a  pressure  force  in 
the  vicinity  of  the  cathode  surface.  The  influence  of  the 
density  gradient  on  the  plasma  particle  motion  is  more 
difficult  to  assess  with  the  given  information.  The  data 
presented  is  static  while  the  effect  of  the  density  gradient 
is  dynamic  in  nature. 


W  «  50  eV,  m  -  1 ,  f  - 1 


t  (nanoseconds) 

Figure  8.  Delay  time  variation. 


W  -  50  eV,  m  ■  1 ,  f «  1 


t  (nanoseconds) 

Figure  9.  Gun  density  variation. 


The  carbon  atoms  forming  the  neutral  cloud  do  not  "feel" 
the  effect  of  the  electric  field.  Therefore,  the  ionization 
of  the  neutral  particles  is  not  affected  by  the  density  < 

gradient.  The  ionized  particles,  on  the  other  hand,  will  be 
accelerated  by  the  electric  field  toward  the  cathode  or  anode 
depending  on  the  charge  and  direction  of  the  field  at  that 
point  in  space. 

The  electric  field  profile  due  to  the  density  gradient 
will  be  of  the  form  presented  in  Figure  10  where  negative  E 
represents  the  direction  toward  the  cathode  and  positive  E  is 
toward  the  anode.  The  data  for  the  electric  field  produced 
by  the  plasma  density  gradient  is  presented  only  between  the 
cathode  surface  and  the  first  annular  region. 


E  (V/m) 


- 

Bulk  Plasma 
Region 


r  (m) 


Figure  10.  Electric  field  profile. 


Figures  11,  12,  13,  14,  and  15  show  that  the  electric 
field  due  to  the  density  gradient  rapidly  saturates.  The 
reason  for  this  is  that  Equation  15  becomes  constant  when  the 
ratio  D  becomes  large  for  a  fixed  electron  energy.  Equation 
15  then  becomes  a  function  of  the  distance  over  which  the 
gradient  is  calculated.  This  is  essentially  a  function  of 
the  expansion  time  of  the  neutral  carbon  cloud  as  shown  in 
Figure  14.  In  addition,  the  energy  of  the  plasma  also  affect 
the  electric  field  of  Equation  15  as  shown  in  Figure  13. 


W  «  50  eV,  m  =  1 


Figure  11.  Density  profile  variation. 


The  purpose  of  this  study  was  to  conduct  some 
preliminary  testing  of  Schwirzke's  Plasma  Sheath  Model  in  an 
effort  to  determine  the  effect  of  the  neutral  carbon 
ionization  process  on  the  electric  field  experienced  at  the 
cathode.  Inherent  in  that  attempt  was  the  determination  that 
ionization  was  indeed  occurring.  In  addition,  the  effect  of 
the  plasma  density  gradient  was  considered. 

It  was  shown  that  even  if  modest  levels  of  surface 
contaminants  were  desorbed,  ionization  was  appreciable.  As  a 
result  of  this  study,  a  baseline  level  of  desorbed  particles 
was  established.  Although  the  results  are  approximations 
based  on  ideal  conditions,  this  information  may  be  helpful 
for  future  investigations. 

It  was  shown  that  variations  in  the  delay  time,  plasma 
density  and  the  desorption  level  provided  the  most  dramatic 
effects  on  the  density  ratio  D.  This,  in  turn,  affects  the 
electric  field  and  the  electron  field  emission  current  drawn 
from  the  cathode  surface. 

The  most  interesting  effect  of  the  neutral  particle 
ionization,  is  the  effect  due  to  the  plasma  density  gradient 
which  was  shown  to  produce  significant  electric  fields  near 
the  cathode.  This  effect  could  prove  to  be  significant  in 
future  studies  for  two  reasons.  First,  the  plasma  density 


gradient  produces  electric  fields  on  either  side  of  the 
ionizing  region  which  are  diametrically  opposed.  Second,  the 
expansion  of  the  ions  is  confined  on  one  side  by  the  cathode 
and  essentially  free  on  the  anode  side. 

The  electric  field  produced  by  the  density  gradient  near 
the  cathode  accelerates  the  ions  toward  the  cathode  while  the 
electrons  are  essentially  trapped  within  the  opposing 
electric  fields  near  the  location  of  E  =  0.  As  the  ions  are 
lost  to  the  cathode  surface,  the  electrons  form  a  negative 
space  charge  layer  which  reduces  the  emission  of  electrons 
from  the  cathode.  This  could  allow  the  current  to  be 
transfered  to  the  load.  At  the  same  time,  the  anode  side  of 
the  density  gradient  produces  an  electric  field  which  opposes 
the  electron  motion  toward  the  anode.  Depending  on  the  rate 
of  expansion  and  ionization,  this  electric  field  could  be 
responsible  for  the  observed  flow  of  ions  near  the  anode.  As 
the  dense  plasma  cloud  continues  to  expand,  the  electric 
field  is  reduced  and  the  density  ratio  approaches  unity 
again. 

The  effects  of  the  anode  pressure  gradient  and  rate  of 
expansion  were  not  included  in  this  study.  However,  the 
dynamic  effects  of  the  density  gradient  and  resultant 
electric  fields  could  prove  important. 
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The  basic  notion  of  the  initial  number  and  number 
density  in  each  of  the  five  concentric  anular  rings  is 
presented  in  the  body  of  Chapter  II  and  Equations  5  and  6. 

The  weighting  factor  is  based  on  a  straight  line  curve 
approximation  where  the  number  of  neutral  particles  remains 
constant . 

The  inner  most  region  is  multiplied  by  the  weighting 
factor  f,  which  ranges  as  0.5  S  f  <  2.  Figure  A1  shows  a 
schematic  of  the  annular  sub-regions  around  the  cathode.  The 
equation  of  the  line  then  depends  on  the  slope  and  the  y- 
intercept .  The  y-intercept  is  determined  by  the  factor  f 
multiplied  by  the  number  of  particles  in  the  innermost 
region.  The  slope  is  determined  by  forcing  the  middle  region 
to  always  have  the  same  number  of  particles  as  was  originally 
allocated. 

This  method  of  weighting  is  somewhat  restrictive 
althougth  for  the  purposes  of  this  study,  the  factor  f  is 
easily  identified.  For  example,  for  f  =  1,  each  region  has 
the  same  number  of  particles  while  f  =  2,  immediately  alerts 
the  observer  that  the  inner  region  is  heavily  weighted  with 
particles . 


Figure  Al.  Annular  sub-regions 
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APPENDIX  B 


program  sheath (input, output) ; 

{The  purpose  of  this  program  is  to  calculate  the  ionization 
of  neutral  particles  which  are  desorbed  from  a  cathode 
surface  in  a  plasma  opening  switch  configuration.  It  also 
calculates  the  electric  field  produced  by  the  density 
gradient.  Calculations  used  are  based  on  fundamental  plasma 
physics  equations  where  no  attempt  has  been  made  to  apply 
kinetic  theory.} 

uses  SANE; 

const 

charge  =  1.6E-19; 
emass  =  9.1E-31; 
imass  *  2.004E-26; 
permi  *  8.85E-12; 
cathradius  =  2.5E-2; 
anoderadius  =  5E-2; 
driftvel  =  7E4; 
tvel  =  330; 
swwidth  =  4E-2; 
swlength  ■  9.5E-2; 
delay  =*  1.5E-6; 
iondensity  =  8E19; 


var 

voltage, ivelocity, evelocity, current, selvolt,neutdensity, 

lamda, sigma, ncur, endcur, catharea, anodearea, eondensity, 

ji,  je, x, ratio, neutradius, 

stepcur, irad, n, length, ef lux, neutpart , 

partsinregion, f, left, slope, parts, lam, ionized: real; 

idens,deltat,dist, ri, ro, pos : extended; 

t,  i,  steps, k, j,m: integer ; 

f ilel : text ; 

chrchar; 

nparticles, iparticles, ndensity, 

|  idensity,ep: array  [1..6]  of  real; 

i  match : boolean; 

I 

I  procedure  setcon;  {set  the  initial  constants  for  the  given 

|  parameters} 

1  begin 

ivelocity : *sqrt (2*charge*voltage/imass) ; 

i 

\ 


evelocity :=sqrt (2*charge*voltage/emass) ; 

neut radius :=cathradius+ (delay-swlength/driftvel) *tvel; 

x:=neutradius-cathradius; 

dist :=x/5; 

catharea : =2  *pi *cathradius *  swwidth; 
anodearea: =2 *pi*anoderadius* swwidth; 
neutpart :  =catharea*n*l .  217E19; 

end; 

procedure  cathodedensity;  {determine  the  current  density  and 

electron  flux} 

begin 

ji :=iondensity*charge*driftvel; 
je:= (current- ( ji* catharea) ) / (anodearea) ; 
eflux:=je/charge*lE-9; if  eflux<=G  then  eflux:=l; 
{part/square  meter  nanosecond} 

end; 

procedure  get info;  {assume  a  straight  line  current  profile} 
begin 

write ('Enter  the  anode-cathode  voltage  '); 
readln (voltage) ; 

write ('Enter  the  number  of  monolayers  to  remove  '); 
readln (n) ; 

write ('Enter  the  starting  current  in  kiloamps  '); 
readln (ncur) ; 
current : =ncur  * 1 E  3 ; 

write ('Enter  the  ending  current  in  kiloamps  '); 
readln (endcur) ; 
endcur : =endcur * 1E3 ; 

write ('Enter  the  number  of  steps  to  take  '); 
readln (steps) ; 

stepcur := (endcur-current) / (float (steps) ) ; 

end; 

procedure  setsigma;  {assign  the  proper  cross  section} 
var 

setvolt : integer; 
begin 

setvolt ;=Round (selvolt) ; 
match :=TRUE; 
case  setvolt  of 

20;  sigma:=lE-20; 

30:  sigma :=1.8E-20; 

40:  sigma:=2.1E-20; 

50:  sigma:=2.3E-20; 

60:  sigma:=2.4E-20; 

70:  sigma: -2 . 3E-20; 

80:  sigma:=2 .2E-20; 

90:  sigma:=2.15E-20; 

100:  sigma:*2 .IE-20; 
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200:  sigma :=1.5E-20; 

300:  sigma:=1.2E-20; 

400:  sigma :=9E-21; 

500:  sigma :=8E-21; 

600:  sigma :=7E-21; 

700:  sigma:=6.5E-21; 

800:  sigma :=6E-21; 

900:  sigma:=5.5E-21; 

1000:  sigma:=5E-21; 
otherwise 

match :=FALSE; 

end; 

end; 

procedure  profile;  {write  the  equation  of  the  line} 
begin 

if  fol  then  begin 

left : =f *partsinregion; 
slope:=2*  <1— f ) *partsinregion/x; 
end; 

if  f=l  then  begin 

left : =part sinregion; 

slope :=0; 

end; 

end; 

procedure  openfiles; 
begin 

rewrite (filel, 'data  set  1'); 

end; 

begin 
openfiles; 
match:=FALSE; 
while  not  match  do  begin 
getinfo; 

selvolt : =voltage; 

setsigma; 

end; 

set con; length : -dist ; 

part  s inregi on : =neu tpart / 1 0 ; 

writeln( 'Enter  the  factor  ' ) ; readln (f ) ; if  f<0  then  halt; 
prof ile;pos:=dist; 

write (filel, 'voltage  =  ', voltage: 11) ; 
write (filel, ' length  =  ',x:ll); 
write (filel, 'monolayer  *  *,n:ll); 
writeln (filel, ' factor  =  ' ,f:ll); 
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for  k:=l  to  5  do  begin  {assign  the  particles  to  the  regions} 
nparticles [k] :=slope* (pos-dist/2) +left; 
iparticles [k] :=0; 
pos :=pos+dist; 
end; 

for  i:=l  to  steps  do  begin 

cathodedensity ; ri : =cathradius ; ro : =cathradius+dist ; 

for  j : =1  to  5  do  begin 

ndensity [ j] :=nparticles [  j  ]  / (pi*swwidth* (sqr (ro) - 
sqr (ri) ) ) ; 

lamda:=l/ (ndensity [j] *sigma) ; 
ionized: = (current /charge- ji) *lE-9* (1-exp ( (- 
length) /lamda) ) ; 

iparticles [ j] :=iparticles [ j] +ionized; 
idensity [ j] :=iparticles [ j] / (pi*swwidth* (sqr (ro) - 
sqr (ri) ) )  ; 

nparticles ( j ] : =npart icles [ j ] -iparticles [ j ] ; 
if  nparticles [j] <=0  then  nparticles [j] :=1; 
ri : =ro; ro : =ro+dist ; 
end; 

if  i  mod  5=0  then  begin  {print  out  every  5  nanoseconds} 

write (filel, i : 5) ;write (filel, '  ' ) ; 

write (filel, voltage: 11) ; 

write (filel, current: 11) ; 

for  m:=l  to  5  do  begin 

write(filel, (idensity [m] +iondensity) /iondensity : 11) ; 
end; 

write (filel ,chr (9) ) ; 

ep [1] :=abs ( (2*selvolt) / (2*iondensity+idensity [1] ) * 
(idensity [1] ) /dist) ; 
write (filel, ep  [1]  :  11) ; 
for  m:=2  to  5  do  begin 

ep[m] :=abs ( (2*selvolt) / (idensity [ml] +idensity [m] + 
(2*iondensity) ) * (idensity [m] -idensity [m-1] ) /dist) ; 
write (filel, ep[m] :11) ; 
end; 

ep[6] :=abs ( (2*selvolt) / (idensity [5] +2*iondensity) *idensity[5] 
/dist) ; 

writeln (filel, ep [6} : 11) ; 

end; 

current : =current+stepcur ; 
end; 

close (filel) ; 

readln; 

end. 
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